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Abstract 

We study integrable vertex models and quantum spin chains with toroidal boundary 
conditions. An interesting class of such boundaries is associated with non-diagonal 
twist matrices. For such models there are no trivial reference states upon which a 
Bethe ansatz calculation can be constructed, in contrast to the well-known case of 
periodic boundary conditions. In this paper we show how the transfer matrix eigenvalue 
expression for the spin-s XXZ chain twisted by the charge-conjugation matrix can in 
fact be obtained. The technique used is the generalization to spin-s of the functional 
relation method based on “pair-propagation through a vertex”. The Bethe ansatz-type 
equations obtained reduce, in the case of lattice size = 1, to those recently found for 
the Hofstadter problem of Bloch electrons on a square lattice in a magnetic held. 


Running title: Spin-s XXZ chain with twists 


1 Introduction 


Vertex models and related quantum spin chains are most often studied with periodic 
boundary conditions |^, |^, |^. If the Boltzmann weights underlying the models are ob¬ 
tained from solutions of the Yang-Baxter equation, or /^-matrices, then commutativity 
of the associated transfer matrices follows and leads to integrability. A large class of 
/^-matrices and therefore integrable vertex models with periodic boundaries can be con¬ 
structed via the machinery of quantum groups ia- Several techniques - coordinate 
Bethe ansatz, algebraic Bethe ansatz, analytic Bethe ansatz, etc. ~ have been developed 
for diagonalizing the corresponding transfer matrices. 

More general toroidal boundary conditions than periodic are in fact integrable; the 
condition being that the “twist matrix” specifying the boundary condition is a “sym¬ 
metry” of the i?-matrix |P . Especially interesting are the cases where the twist matrix 
is non-diagonal. In spite of the existence of commuting transfer matrices for such cases, 
their diagonalization by the above mentioned Bethe ansatz techniques seems to fail. 
This is due to the non-existence in such cases of a trivial reference state. 

In this paper we show that the spin-s XXZ chain twisted by the (non-diagonal) 
charge-conjugation matrix can in fact be diagonalized. The technique used is a func¬ 
tional relation method (also known as the “T-Q” method or the “method of com¬ 
muting transfer matrices”) based on the “pair-propagation through a vertex” property 
of the relevant i?-matrices. This generalizes to spin-s a recent result for the spin-f chain 
(or six-vertex model) |0]. 

One way [^, to study the spin-s XXZ chain with periodic boundaries is to hrst 
express the transfer matrix (u) in terms of a much simpler auxiliary transfer matrix 
t^^\u) using the fusion procedure |^. This auxiliary transfer matrix is then diagonal¬ 
ized with the algebraic Bethe ansatz to give the eigenvalue expression and associated 
Bethe ansatz equations for the original spin-s chain. We show that in the presence of 
the twist, there is an analogous procedure. The corresponding auxiliary transfer ma¬ 
trix is then diagonalized with the abovementioned functional relation method. As a 
consquence we obtain the transfer matrix eigenspectrum for the spin-s XXZ chain with 
the charge-conjugation twist. 

The paper is organized as follows: Section 2 is introductory in nature. We describe 
the general theory of integrable toroidal boundary conditions here and make a survey of 
available results for the case of the spin-f XXZ Heisenberg chain. In Section 3 we study 
the integrable spin-s generalization of the XXZ chain twisted by the charge conjugation 
matrix. In particular we relate the associated transfer matrix t^ 2 s‘\u) to the auxiliary 
transfer matrix which is dehned in terms of the i?-matrix Section 4 

is the body of the paper. In Section 4.1 we show that the property of “pair-propagation 
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through a vertex” for the i?-matrix (u) and related spin-^ XXZ chain |I| generalizes 
to arbitrary spin-s. This allows the T — Q relations for this auxiliary transfer matrix 
to be obtained in Section 4.2. For this purpose we required a conjectured property for 
two families of matrices Ql{u) and Qr{u), which is partially proved in Appendix A. In 
Section 4.3 we obtain the eigenvalue expression and Bethe ansatz-type relations for the 
auxiliary transfer matrix, leading to the corresponding expressions for the spin-s XXZ 
chain. In Section 5 we discuss the contents of this paper in a general context. 

After this work was completed, we became aware of an interesting connection be¬ 
tween integrable non-diagonal twists for the auxiliary transfer matrix in the 

case of lattice size N = 1 and the so-called Hofstadter problem of Bloch electrons on a 
square lattice in a magnetic held In Appendix B we make contact with the latter 
work; in particular, recovering their Bethe ansatz results derived using a completely 
different method. 


2 Integrable models with toroidal boundary condi¬ 
tions 

Let us recall a few fundamental facts concerning the quantum inverse scattering method 
and the construction of integrable models The principal ingredient is the Yang- 

Baxter equation 


Rl2{u)Ri3{u + v)R23{v) = R23{v)Ri3{u v)Ri2{u) ( 2 . 1 ) 

which is an equation in the space Vi®V 2 ® V 3 , with y = Y ~ (D”, R{u) E End(lZ ® V) 
and Ri 2 {u) = i?(M)®l, R 2 z{u) = l®R{u), etc. Given an i?-matrix, which by dehnition 
is a solution of 0) one constructs a monodromy matrix 

T (n) = Rai{u)Ra2{u) ■ --RaNiu), (2.2) 

which is to be thought of as an operator-valued matrix in the auxiliary space 14. A 
consequence of (|2.1|) is the intertwining relation 

Ri 2 iu -v)T (u) T (v) =T (v) T {u)Ri 2 {u - v) (2.3) 

for the monodromy matrices. The indices 1 and 2 here label different auxiliary spaces 
and are not to be confused with those appearing in the dehnition (p.2|) which label 
different quantum spaces. The important object 

tp{u) = tr„ f (n), (2.4) 


2 





has the interpretation of a row-to-row transfer matrix for an n-state vertex model with 
periodic bonndary conditions and forms a commuting family. 


We consider an extension of (|2.4|) to more general toroidal boundaries; namely by 
studying the transfer matrix 

t{u) = tia f (m) F (2.5) 


where F is a square matrix of dimension diml4. The matrix F is called a “twist matrix” 
and the boundary conditions twisted. The following simple result shows that for certain 
twist matrices F commutativity of transfer matrices is retained: 


Proposition 1 Let the twist matrix F he a “symmetry” of the R-matrix, i.e it obeys 
[R{u),F ^ F] = 0. Then the transfer matrix ^2.^ ) satisfies [t{u),t{v)] = 0 and we say 
F is integrable. 


Proof From the dehnition, we have t{u) t{v) = tri 2 FT {u) T (n) F- Using the 

/iTi-ra ^ ^ 

interwining relation (|2.5|) to change the order of T {u) T (n) we obtain 


1122 1 21 2 

tri 2 FT {u) T (v) F = tri 2 F Rfi{u - v) T {v) T {u)Ri 2 {u - v) F 

= tii 2 Rff{u -v)T {v) T {u)Ri2{u - v)F (g) F 
2 1 

(symmetry) = tri 2 F ;^2 ("^ — v) T {v) T {u)F ® FRi 2 {u — v) 
= tri 2 T (n) T {u) FF = t{v) t{u) 


A quantum spin chain can be associated with the transfer matrix (|2.5|) in the usual 
way. For simplicity we assume the F-matrices are regular, with i?(0) = cP, where P 
is the exchange operator with matrix elements = 6^6 ^. Given the multi-indices 
(oi,..., Oat) and {Pi,, Pjsi) we have from (p.5|) 


,, xPn TP01 

Similarly, assuming F is invertible, we have 

C---€:z\(F-vmLR') 

From this we conclude that the quantum spin chain has Hamiltonian 

1 \ 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 
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with 


1 N 

hi,j = -Rij(0)'Pii, =F~^ hN,i F, (2.9) 

P 

where F acts non-trivially only in the p-th slot of l/^ ® ® V^. When the twist matrix 

F is the identity, the spin chain with periodic b.c.’s is recovered. 

Example: Spin-f XXZ chain 

The six-vertex model i?-matrix is given by 

/ sinh(M -|- A) \ 

^ sinhw sinh(A) 

^ ' sinh(A) sinliM 

\ sinh(M -|- A) / 


and admits the following two types of integrable twist matrices: 


(i) F 


a 0 \ 

op)' 


(ii) F 


0 a \ 
P 0 


where a, P are completely arbitrary. When a, (3 are non-zerothe corresponding spin-^ 
XXZ Hamiltonian is given by 


1 ^ 

^ _ _ \ ^ 

sinh(A) ^ 

with boundary conditions 

(i) ^7V+l 


^j+i 




cosh(A) 




p + 

a 

= 

a 


P - 

a 

= 1 
a 



1 , cr^+i = al 


( 2 . 10 ) 


( 2 . 11 ) 


For boundary conditions of type (i), the total spin operator P2j=i commutes with 
the Hamiltonian; the trivial reference state is available and the corresponding transfer 
tOpu) matrix can be diagonalized with the algebraic Bethe ansatz [^, with very little 
difference from the periodic case. The transfer matrix eigenvalue expression is given by 


A0 \u) = asinh(M-I-A)^ J)[ 


sinh(M — Uj — A) 


sinh(M — Uj) 

n ■ ^ . M PA siuhiu - Uj + X) 

-F Psmh{u)^Y[^\- - 

smh(-u — Uj) 


*When one of a,/3 is allowed to be zero, the vertex model has fixed boundaries; it does 
not seem possible to obtain a local Hamiltonian in this case. Furthermore, the eigenvalue 
expression for such type (ii) twists is unknown. 


( 2 . 12 ) 
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in the sector with total spin N — 2n, with the Bethe ansatz equations for Uk determined 
from analyticity of 

On the other hand, boundary conditions of type (ii) leave only a Z 2 -invariance in 
the eigenspectrum. For the case a = /? = l[| (which can be considered as anti-periodic 
boundary conditions on the vertex model), the problem was recently solved using 
the functional equation method of Baxter |Q. The transfer matrix t{u) has eigenvalue 
expression 


h(u) 


smh(«+Af 

/=i smh[i(M-Mj] 

si„h(«)» n + A)] 

,Vi smh[yti - Uj] 


(2.13) 


General boundary conditions of type (ii) are in fact very closely related. This is due to 
the relation 



(2.14) 


between the characteristic polynomials for their corresponding transfer matrices, leading 
to the result 



(2.15) 


Type (i) boundary conditions with aj3 = 1 were studied in Ref. [^]: At criticality 
the model gives rise to a conformal held theory with central charge c < 1, with the 
twist parameter determining the value of c. All the above toroidal boundary conditions 
have in fact been studied in Ref. fl^ (by direct numerical diagonalization in the case 
of type (ii) twists), with conjectures given for their complete operator content. In their 
language, the type (ii) boundary conditions share the same physical properties because 
they are in the same conjugacy class of 0(2), being the global symmetry of the inhnite 
chain. 


3 Spin-s XXZ quantum chains 


The integrable spin-s generalization of the XXZ spin chain and related R-matrices have 
been studied by many authors; amongst them o, 0, l, 0 0 ig 0 0 


In 


particular, the spin-1 Hamiltonian was hrst constructed in Ref. |^. Bethe ansatz 


equations for the spin-s models with periodic boundary conditions were obtained and 
studied in Refs. B0S. For integrable boundary conditions associated with diagonal 
twists, Bethe ansatz equations were obtained and studied numerically in Ref. [|^ and 
'I'We refer to this as the charge-conjugation twist 
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analytically for spin-1 in Ref. [|T9|. General toroidal bonndary conditions for the spin-1 
case was also investigated by direct nnmerical diagonalization in Ref. fpaper, 

we will mainly follow the approach of Ref. to constrnct the integrable spin-s chain. 

The transfer matrix for the spin-s XXZ chain is associated with the R-matrix 

^(2.,2.) 

(u) which acts in a tensor prodnct of two spin-s representations of l/g(su(2)). 
We will stndy it via its constrnction by the fnsion procednre. For this pnrpose we will 
reqnire the R-matrix acting in the tensor prodnct of the spin-Z/2 and spin-m/2 

representations of Uq{su{2)) for arbitrary l,m. These can all be constrncted iteratively 
from R*'*’^^(m), whose explicit form is given by 

= sinh -I--I- ® -I- sinh(A) S~^ a~ + S~ ^ . (3.1) 

Here and 5^ are generators for Uq{su{2))\ 

= ±25^, 

In the spin-//2 representation these operators take the forms = diag(/, I — 2,, —1) 
and 

/O /(I) \ 


= [s-y 


sinh A 


V 


/(O 

0 


where f{j) = Y^sinh(jA) sinh[A(/ -|- 1 — j)]- The operators and cr^ are in the spin-1/2 
representation. 

The R-matrix degenerates at n = A and u = —/A: 

R('’b(A) = R('>^)(-/A) = (3.3) 


where and are matrices and are projection operators onto snbspaces of 

dimension I -|- 2 and I, respectively. Using this property and the fnsion procednre 0 it 
is possible to establish the key relations [§ 

+ A)R&™V)^f2 = 

RfuT\^) 0 

* sinhMsinh[M + {m + 1)A]R<]^2^>™^('^ + 2A) 

/ sinh(M-|-mA)R|^j^^) 3 ™^(-u) 0 

y * sinhMsinh[n + [m + 1)A]R<]^2^>™^('*^ + 2A) 

/ sinh(M-|-mA)R| 32 ) 3 ™^(-u) o\ 
y * sinhMR<y2^y™^(n-|-2A) j 


^ if / < m, 

I if Z = m, (3.4) 
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Expressions for the matrices and are given in [|^. The spaces V(i 2 ) 

and E<i 2 > in eqnation ( |3.4|) are the symmetrized and anti-symmetrized components 
respectively of Vi ® V 2 . The relations (|3.1|) and (|3.4|) allow the i?-matrix (u), 
throngh which we obtain the spin-//2 XXZ chain, to be obtained recnrsively. A non- 
recnrsive dehnition is also available [§. 

Thus far, all our considerations have been purely local in character; in particular, 
they are independent of boundary considerations. We will now specialize to the bound¬ 
ary conditions specihed by the non-diagonal twist matrix (the charge-conjugation twist) 
^(0 = , which is the (/ -|- 1) x (/ -|- 1) matrix with 1 along the anti-diagonal and 

0 elsewhere. It has the properties 

^(0^±^(0 = s^, = -FT (3.5) 


Using this it is easy to see from ( p.l|) that ® = 0. The recursion 

relations (^) then imply that F^^^ is a symmetry (in the sense of Proposition 1) of 

Dehne the monodromy matrices 


r<'>(u) = fli'yV) 


p{t2s) / \ 

■ ■ ■ ^aN (“) 


( 3 , 6 ) 


with a spin-1 auxiliary space and a quantum space which is an Wfold tensor product 
of spin-s representations, and transfer matrices 

= tTaT^^\u) F^^\ ( 3 . 7 ) 


By Proposition 1 (or rather, the obvious generalization to the case of non-isomorphic 
quantum spaces), these form a commuting family: {v)] = 0 for all l,n. 

With respect to the same decomposition as in (|3.4[) we also have 




F(^+i) 0 \ 

0 _F(^-b ) • 


(3.8) 


Consider now the product From the dehnition (|3.7|) and commuting 

the F-matrices appropriately this becomes 


tr.t RiF(u + ■' ■ATP + 


which can be written in the form 

tr., d')FW Di'Fj (AAi7'’('“+ A)dU’('“)r>2^") 
• ■ • (oi'MTF + A)FSt'*>(«)Fr‘) ■ 

By application of eqs (|3.4|) and (p.8|) we arrive at the following result; 
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Proposition 2 The transfer matrices t?^\u) obey the recursion relations 
+ X)t^t\u) = 

f t\‘^i{u) — sinliM^sinh[M + A(2s + l)]^t\‘^^{u + 2A) I < 2s 

< sinh(M + 2sA)^4+*i (m) — sinh sinh[M + A(2s + {u + 2A) I = 2s (3.9) 

[ sinh(M + 2sX)^t^^}{u) — sinh. + 2A) I > 2s 

Define the shift operator by z~^f{u)z = f{u + A) for all functions f{u) and the 
function S^^u) by S'^^u) = sinhM'^sinh[M + A(m + l)]^. Then the generating function 
for the transfer matrices t^j^^\u) is given by 

1 — ztf^\u) — z^S‘^^\u) = 

2s oo k—1 

flz¥k'\u)+ Y. n sinh(M +A/)^4 ^'*^(m). (3.10) 

k=0 l=2s-\-l l=2s 


(on the 


This can be proven by multiplying both sides by ^ — zt^^\u) — z‘^S‘^^\u) 
right) and using the recursion relations ( p.9|) to show that all coefficients of z^ for k>l 
vanish. The corresponding generating function P] for periodic b.c.’s differ from (|3.1(]|) 
in only the sign of the coefficient of the z'^ term on the Ihs. 

The transfer matrix t 2 s^\u), out of which the Hamiltonian for the spin-s XXZ chain is 
constructed, can be obtained in terms of the transfer matrix t^i^\u) either using the re¬ 
cursion relations ( p.9|) or the generating function (|3.10|) . It is this transfer matrix which 
is manageable; for instance, because the auxiliary space for t^‘^^\u) is two-dimensional, 
the functional relation method of Baxter []^ can be generalized in a straightforward 
manner to obtain the eigenvalue expression and associated Bethe ansatz-type equa¬ 
tions. 

At M = —(2s — 1)A the i?-matrix becomes proportional to the exchange 

operator P, with proportionality constant Ojlisinh(jA). According to Eq. (|2.8|) the 
Hamiltonian for the spin chain is then given by 


du 


d 


12^)/ 


n=—(2s—1)A 


2s 


nsinh(jA) ^ 

i=i 




(3.11) 


with two-site interaction = R^^^'‘^^\—{2s — l)X)Pij. The Hamiltonian in terms of 

Uq{su{2)) generators is very complicated, but using Eq. ( |3.5|) the boundary conditions 
can be simply written as 


QZ _ QZ ci _ CT 


(3.12) 
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4 Transfer matrix diagonalization 


In this section we will obtain the eigenvalue expression for the transfer matrix 

from which the eigenvalue expressions for the transfer matrix and the 
Hamiltonian ( |3.11| ) follow. For this purpose we will require explicit expressions for the 
non-zero matrix elements of which can be obtained by interchanging S'* and 

(T* in (^.11): 


where 


= ay(ti), 

/?“•"*'(«)? hi = hW, 

= hW. ( 4 . 1 ) 


aj{u) = sinh[M -|- (2s -f 1 — j)A], 

bj{X) = ^sinh(jA) sinh[A(2s -|- 1 — j)]. (4.2) 


4.1 Pair-propagation through a vertex 

A key ingredient of Baxter’s “method of commuting transfer matrices” [|l| is the property 
of “pair-propagation through a vertex” for the i?-matrix; ie. the existence of vectors 
n™ and n”*** such that:0 

( 4 . 3 ) 

for a, G {1,..., 2s -1- 1} and p, G {1, 2}. Here n™ are entries of the vector n™ etc. 
For the i?-matrix the set of equations ( [4.3|) can be written as 

aj{u)vf h^ + bj_^{X)vf_^hf = (4.4) 

a2s+2-Au)vfhi^ + b,{X)vf_,,h^ = (4.5) 

where j runs from 1 to 2s -|- 1. The first thing to note about the system of equations 
is that it is linear in and Therefore, consistency of ( |4.4| ) and (|4.5|) 

formulated here for i?-matrices acting in 


(O), (T^) 
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requires the vanishing of the (4s + 2) x (4s + 2) determinant 


aihfi 

h2sK 


a2s+iK biK 


h2sK 





An examination of small values of s leads us to the following conjecture: 


(4.6) 


Claim 1 Define 


7 in,out /1 in,out 
^2 / ^1 


Then the determinant ( U-(^ vanishes if 


where a = —s, —s + 1,..., s. 


^OUt 

,^in 


= e""”, 


(4.7) 


We do not have a proof for general s, but this result has been important in leading 
to Proposition 3 below, which is a stronger result and will be proved. For that purpose, 
we will require a few dehnitions: Let q = and dehne the g-numbers, factorials and 
binomial coefficients, respectively, by 



m 

n 


qn _ q-n sinh(nA) 

q — q~^ sinh(A) 

[n] [n — 1] ■ • • [1], with [0]! = 1 

[m ]! , „ 

(m, n G Z>o, m>n). 


The g-binomial coefficients satisfy the following relation (a g-analogue of the binomial 
theorem) 


2s+l 


k=0 


2s + 1 
k 


= n (l + ag^^) , 

j=-s 


(4.8) 


which will be used later. Introduce the functions f^'P\a]u) = fjf\o",u), with 



i-e 


2Xa\k 


p 

k 


k p 

J([ sinh[u + (j — 1)A] sinh[M + (2s + j — p)A]. (4.9) 

j=l j=k+l 


They have the useful properties summarised in the following: 
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Lemma 1 Let ajiy) and bj{X) he defined as in ISD- For any a E €, the functions 
f^\a;u) defined in iM) satisfy 

f^P\a;u) + e‘^^'^a2s+2-piu)f^^~^\a;u) = ai{u)f^P~^\a;u-X), (4.10) 

ap+i{u)f^P\a-,u) + e‘^^'^bp{ufif^P~^\a;u) = f^P\a-,u - X)ai{u), (4.11) 

and 


ai{u)a2s+i{u)f^^ ^\a]u-X)f^^ ^\a]u + X) 

= aj{u)a2s+2-j{u)f‘'^~^\a-,uf - bj^fiX)^f^^\a-,u)f‘'^~^\a-,u). (4.12) 

If a takes values only in the range {—s, —s + 1,..., s}, then we also have 

f^^^\a; u-X)= u). (4.13) 

Proof These relations can be proved using the dehnition of f^\a;u). For instance, 
Eq. (^4.1U|) boils down to verifying 


p — 1 

k 


sinh(M — A) + 


p — 1 

k-1 


sinh[M + (p — 1)A] = 


p 

k 


sinh[M + {k — 1)A]. 


The proof of Eq. ([4.13|) is more interesting; the difference between the two sides can be 
shown to be 


2s 2s+l 

n + U - 1)A1 i: (-e^")" 

j=l k=0 


2s + 1 
k 


Upon application of the g-binomial theorem (|4.8|), the sum above simplifies to 


S 

n(i- 


,2A((T+j) 


j=-s 


from which the result immediately follows. 

We are now ready to state and prove the main result in this subsection: 


Proposition 3 Let a E {— s, — s + 1, 

relations 


, s} and f^\a;u) be defined as in l U-^) - The 


F' 


,out 


V 


in 

i+i 






ai{u) 

1 

W) 


lifivf /(i-b(a;^-A) 
hf^^ /b“i)((T;M)’ 


(4.14) 

(4.15) 

(4.16) 


with arbitrary hf^f, hfi and solve the equations MU and M)- In other words, 

the R-matrix does indeed have the property of pair-propagation through a 

vertex. 
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Proof The proof is by substitution of ( 4.14|) , ( 4.15 ) and ( 4.16|) directly into (^) and 
(^4.5|) . For “generic” values of j, these equations reduce to ( [4.10|) and ([4.11|) which we 
have indicated above are true. For the “end values” j = I and j = 2s + 1 equations 
(^4.4|) and ( |4.5|) , respectively, degenerate by virtue of the vanishing of bo^X) and & 2 s+i(A). 
The former case can be trivially seen to be true; the latter is true by virtue of ( [4.13| ). 
For the special cases of a = ±s, the function f^P\a;u) simplihes to 


sinh[(2s — j + 1)A], 


(4.17) 


with corresponding simplihcations in the expressions for n™ and Vj in Proposition 3. 
For the spin-4 case, we recover the six-vertex model results of by substituting Eq. 
([4.17|) into Proposition 3. 


Let us indicate explicitly the u- and a-dependence of the vectors v^^{u) and n 


out / 


u 


Equation ( 4.16 ) is a recursion relation for the elements v™{u)j^ given the initial value 
v'^{u)i. If we choose this initial value to be unity, then it is evident that 


\u - X)j = n“(n). 


^)(cr;M - A) 


/L i)(cr;M) 

It can then be shown from ([4.15|) that is simplify related to by 


h'm 


[u 


A). 


(4.18) 


Proposition 3 can be interpreted in the following way: It guarantees the existence 
of matrices 


P = 




P' = 




such that the generalized similarity transformation 


(P® l)-^M2’""^(n)(P'0l) = 




a[u) 

7(n) 


(4.19) 


f3{u) 

5{u) 


(4.20) 


results in a matrix which is upper block-triangular when acting on a certain vector in 
to be specified. To see this statement, pre-multiply ([4.2U|) by P ® 1 on both sides 


and act on 1 ® v^^{u). We end up with a 2 x 2 matrix equation with entries in 
Component-wise, it reads 


R{u)\h^^{u)k P(M)iJn“(n)fe 
R{u)l]P^{u)k P(M)gn”(M)fc 

/ t, out . \ , 



* 


a{uy^{u)j (5{u)P^{u)j 


(4.21) 
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A comparison with (O) immediately indicates that 

7(m)C(«) = 0. 


(4.22) 


The action of a{u) and 5{u) on the vector whose components are determined 

from the recursion relation (|4.16|) with u“(m)i = 1, is determined by examination of 
the explicit forms for u“(m) and given in Proposition 3: With the help of (|4.18|) 

we hnd that 


a{uy^{u) = ai{u) - A), 


(4.23) 


whereas by taking determinants on both sides of ( [4.21|) , together with the result (|4.12|) , 
we hnd that 

hrMet(P' 




a2s+i{u) u“(m + A). 


(4.24) 


h'l det(P) 

All the considerations in this subsection are local in nature. The results are directly 
applicable to the transfer matrix with periodic boundaries, generalizing the spin-A (six- 
vertex model) results of |I| and providing an alternative to the algebraic and analytic 
Bethe ansatz p|. In the next subsection we will specialize to the boundary conditions 
specihed by the twist matrices P^^^ dehned in Section 3. 


4.2 Diagonalization of the auxiliary transfer matrix 

The main idea of the technique to diagonalize the transfer matrix {u) is to use Eqs. 
(^4.20|) and (|4.22|) - (^4.24|) to construct functional equations satished by We 

note that an attempt to do a generalized algebraic Bethe ansatz calculation, along the 
lines of for the eight-vertex model, should also probably use Eqs. ( |4.2U|) and ( |4.22|) - 
(^4.24| ) as a starting point. 

Let us rewrite the auxiliary transfer matrix in the form 

= tr P, . (4.25) 

where the matrices Pk are understood to act non-trivially only in the auxiliary space 
Vn. These matrices are chosen to be 


Pv+i — 
P2={ 


/ 1 * 
I r * 

1 


Pn = 


re 


1 * 
—2A(Tjv ^ 


N 

=2' 


* 

* 


Pl = 




re 


(4.26) 


where cxj, for i = 1,..., is any set of variables taking values in {—s, —s -|- 1,..., s}. 
If, in addition, r is chosen to satisfy 

/ N \ 


= exp 2A ^ (Tj 


(4.27) 


2 = 1 
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then, by direct verification, we have the relation 


P, 


-1 

Af+l 


0 1 
1 0 


Pl = 


1 

0 


* 

det(Pi) 


(4.28) 


det(Piv+i) 

Let us denote cr = (ui,..., ajsi). From the results ([4.19|) to 


we deduce that 


tf'\u) Vaiu) = -ai{u)^ya{u - A) - ra 2 s+i{u)^ya{u + A), 
r 

(g)iV 


(4.29) 


where ya{u) = v^^{u) ® ■ ■ ■ ®v^{u) is a vector in , with the components of 


v^J{u) determined from the recursion relations 

= 1, 


.,{k) 

o-fe 


{u)j+i = 


2A((tih ho-fc) f (^)(^ ak ', u ) 




(4.30) 


Note that the r’s in Eq. ( |4.29| ) are cr-dependent because of Eq. (|4.27|) ; they can be 
removed by re-scaling ya{u) by a n-dependent factor. The result can be summarised 
thus: 

Proposition 4 Let be the transfer matrix defined by (|J. 1[). The vector 


qaiu) = 


(4-31) 


where v^^fu) is determined by 4^. 5*61 ) with r = exp(A <^0? satisfies the equation 

qcr{u) = ai(n)^gcr(M - A) - a 2 s+i{u)^qa{u + A). (4.32) 


The vector qaiu) in Proposition 4 is well-defined for any choice of cr = (cxi,..., ctat), 
with (jj G {—s, —s -|- 1,. .., s}. Therefore the matrix Qr{u), whose columns are formed 
from the collection of such vectors ((2s -|- 1)^ of them altogether), satisfies the relation 

tf"\u) Qr{u) = ai{uYQR{u - A) - a2s+i{uYQ r{u + A). (4.33) 

We will now require the following property of (u) which is due to crossing-symmetry 
of the corresponding i?-matrix: 

Lemma 2 The transfer matrix {u) satisfies the relation 

- 2As) = (4.34) 
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Proof The i?-matrix (u) has the crossing relation 

=V - 2As) P, (4.35) 

with V = a~^ — a~, which can be verihed directly from Eq. (0). By the properties of 


the trace, we can write the Ihs of (^4.34|) as 

tr,*“*--*- - 2As) ■ ■■Ri]^^^\-u - 2As) 

= tr„ - 2As) • ■ - 2As)) , 


which upon application of Eq. (|4.35| ) becomes 

= ujv-'Ri\-^‘>(u)V- 


r-1 d { 1 , 2 s )^^,^ \ I y-i ^-1 ^( 1 ) _ 


Upon using the properties (U~^)^ = —1 and = F^^\ the assertion follows. 

With the relation (|4.34|) , together with ai{—u — 2Xs) = —a 2 s+i{u) and a 2 s+i{—u — 
2As) = —ai{u) which follow from the dehnition ([4.2|), we conclude from Eq. (|4.33|) that 


Ql{u) = ai{u)^QL{u - A) - a2s+i{u)^Q l{u + A), 


(4.36) 


where we have dehned Ql{u) = ^Qr{—u — 2As). The matrices Ql{u) and Qr{u) obey 
the crucial “commutation relations” 

Ql{u)Qr{v) = Ql{v)Qr{u). (4.37) 

This follows from the following result which we formulate as a conjecture: 

Conjecture 1 For any choice of ai, a[ taking values in {—s,..., s}, the inner product 

^qa'{-u-2Xs) ■ qa{v) (4.38) 

with the vector qaiu) defined in Proposition 4, is symmetric in iu^v). 

Due to the tensor product structure of the vector qcr{u), the inner product (^ 
assumes the form 

fiafi-u - 2Xs) ■ qa{v) = i-D-E"ifj 


k=l 




(4.39) 
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where 


G^^)^^{u,v) = - 2\s) ■ v^^^iy). 


However, from the recursion relations ([4.30|) together with r = we have 




61(A) ■■■6,(A)’ 

where denotes Xil^+i giving rise to 

^ g2jAX''(0-+cr') 


(4.40) 


G (m, v) = ^ 


^o(6i(A)---6,(A))^ 


-u - 2sA). (4.41) 


If we restrict crj,cr' to take values only in the set |±s| then (n, n) in Eq. (|4.39|) 
simplihes and it can be proved (see Appendix A) that the inner product (|4.39|) is indeed 
symmetric in {u,v). For s = | this is, of course, the whole story [|^. For general s we 
have not been able to construct a proof of Conjecture 1 but we have conhrmed the 
result on a computer for a range of values of s and N. 

We now require a technical assumption; 


Claim 2 The matrix Qr{u) is non-singular for generic values of u. In particular, it 
can be inverted at the point u = Uq. 


Dehne now the matrix Q{u) = QR{u)Qf^{ uq). With the property (|4.37|) and the above 
technical assumption (which is standard in the present method |^, but is presumably 
difficult to prove), we arrive at the key result that 


t^i’'\u) Q{u) = Q{u)t^i^\u) 


= ai{u)^Q{u - - a 2 s+i{uy'^Q{u + X), 


N/ 


(4.42) 


together with Q(u)Q(v) = Q(v)Q(u). In other words, we have constructed a matrix 
Q(u) which commutes with the transfer matrix ty^\u), and with Q(v) for v different 
from u. This allows us to choose the diagonal representation for all the matrices appear¬ 
ing in the functional equation ([4.42| ), which then becomes a relation for the eigenvalues 
of the transfer matrix. 


4.3 Bethe ansatz-type equations 

We can deduce the functional form of Q(u) in its diagonal representation (written also 
as Q{u), with abuse of notation) by examining Eq. (|4.42|) in the limits u —^ ±cx). 
From Eq. (|3.1|) the dominant term in the large u limit in each entry of R^’'X(^u) is ~ e“. 
Therefore from Eq. (p.7D the dominant term in each entry of ty^\u) is ~ (rather 
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than 6^“ because of the twist matrix On general grounds the large u behaviour 

of Q{u) is argued to be ~ Subsituting into Eq. ([4.42|) we find that r must be sN. 
A similar argument holds for u —»• — oo and we are led to the form 

Q{u) ~ + • • • + 


which can be parametrized alternatively as 


2sN 

Q{u) = n sinh 
i=i 


u — u 


(4.43) 


after taking out an irrelevant overall u-independent factor. The eigenvalue expression 
for the transfer matrix can then be written as 

Ap'(u) = + 

,=i smh[i(n-nj)] 


- sinh(n)^Tf 


(4.44) 


The requirement of analyticity of A^^*^(m) then imposes the following conditions on the 


‘Bethe ansatz roots” Uj: 


( sinh(Mfc + 2sA) ^ ^ sinh[4 (m^ - Uj + A)] ^ ^ 


y sinh(-UA;) 


sinh[4(Mfc - Uj - A)] 


(4.45) 


In the spin-4 case, we recover from (f4.44|) and (^4.45|) the eigenvalue expression and Bethe 
ansatz equations obtained in 0. For general s, substitution of (|4.44| ) into the generating 
function ( p.lOj ) yields the eigenvalue expression for the transfer matrix [u) in terms 
of the Bethe ansatz roots Uj, namely 


i-i 


Af'w ^ Muf n X 


2sN 

n 


where 


sinh[^(-u — Uk + /A)] sinh[^(M — Uk — A)] 
sinh[4(M -Uk+ jA)] sinh[4(M - Uk + {j - 1)A)] 

( T I h TT 


I-I i-1 

= n sinh(M -|- pA) sinh[M -|- (2s -|- p)A]. 

p=j p=0 


+ 


(4.46) 


(4.47) 
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The eigenvalue expression for the spin-s XXZ Hamiltonian ( 3.11 ) follows from 


A 25 ^^(m) in Eq. (^4.46|) : We find the following expression 


A. 

du 


logA2'^(M) 


w=—(2s—1)A 


2s 2sN 

NY,cothU^) + T, 

j=l k=l 


sinh(sA) 

sinh[|(2sA + Uk)] sinh(AMfc) ’ 


(4.48) 


with the Bethe ansatz roots Uk determined from Eq. (|4.45|) . 

We close this section with a few remarks on more general non-diagonal twists for 
the transfer matrix. Namely, consider 


a 


ftHu) 



with 

f(^) = 1 

0 a \ 

J 0 )■ 

We find the relationship 




det (m) — oc det {tA"’ ~ Aj 


(4.49) 


(4.50) 


between the characteristic polynomials for the two transfer matrices, giving the eigen¬ 
value expression 


A)' 


(2s) 



(4.51) 


with A^^^^(-u) given by Eq. ( [4.44| ). The twist matrix generates recursively the 
twist matrices with the matrix elements of given by f/j = 
through the decomposition 


which replaces Eq. ( 3.8[) when a, (3 are different from 1. The auxiliary transfer matrix 
is then relevant to the twist matrix for the spin-s XXZ chain corresponding 
to boundary conditions 

3 Oi 

^N+l — ^N+l — 1 ^N+l — (4.53) 

generalizing Eq. (p.llD to arbitrary spin-s. By Eq. (|4.51|) the Hamiltonians correspond¬ 
ing to these boundary conditions share the eigenvalue expression Eq. (3.481) and Bethe 
ansatz equations Eq. (3.45|) . 
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5 Discussion 


In this paper we have studied integrable toroidal boundary conditions for vertex models 
and related spin chains, concentrating on boundary conditions specihed by non-diagonal 
twists. In particular, we have applied the “T-Q” functional relation method to obtain 
the transfer matrix eigenvalue expression and Bethe ansatz-type equations for the spin- 
s XXZ Heisenberg chain, thereby generalizing a recent result |]^ for spin-f. A detailed 
analysis of these Bethe ansatz equations should conhrm various conjectures available 
[^, 1^ for the critical properties of these models. 

The primary reason why the abovementioned functional relation method is appli- 
cable is because the i?-matrix giving rise to the transfer matrix t\ has the “pair- 
propagation through a vertex” property. Such considerations make sense in the hrst 
place because of the two-dimensional nature of the auxiliary space 14 associated with 
this transfer matrix. For more general /^-matrices such as those given in Refs. ii, 
where the auxiliary space is no longer two-dimensional, the appropriate generalization 
is not so clear. Therefore the diagonalization of the transfer matrices for such models 
in the case of non-diagonal twists remains an open problem. 
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Appendices 

A Partial proof of Conjecture 1 


In this appendix we study the inner product *qcr'{—u — 2Xs) ■ qcr{v) of Conjecture 1, for 
the case when Uj, are allowed to take only the values ±s. In this case, the function 
simplihes according to (|4.17| ), accompanied by a simplihcation of 
in ( |4.41| ) to 


{u,v' 


2s 

= E 

j=0 


2s 

j 


^jXj:\a+CT') .(2<A-^+^) 


(A.l) 
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Here 


2s 

J 


is the g-binomial coefficient defined in the paragraph preceding Lemma 1. 


By an application of the g-binomial theorem ( [4.8|) we obtain 


s-l/2 / 

Ga^'aiu,v)= n 1 + e 

j=-s+l/2 \ 


2j\+XJ2\a+cr') 2 <t'^\- 


(TuV , 


(A.2) 


and therefore 


^qa'{-u - 2As) ■ qa{v) = 


N 


s-l/2 

n n 

k=l j=-s+l/2 




(A.3) 


Snppose now there exists a pair cTp, snch that ap + a'^ = 0. Then the terms in Eq. 
(|A.3|) which involve cXp can be seen to be symmetric in (M,n), whereas the terms which 
do not involve cXp can be seen (after some relabelling) to be expressible in the form (|A.3|) 
with N N — 1. Hence we need only consider the case where Uj = a[ for all i, which 
simplifies to 


Gn{(Ti, ...,ajv)= *qcr(-u - 2As) ■ qcr(v) = e 
k=lj=-s+l/2 V 


E n 
i = . 


= k + l 


E k — 1 
i=l ^ 


(A.4) 


If (jfe + CTfc+i = 0 then from Eq. (|A.4|) we can write Gn{o'i, ... ,ak, —cxfc,..., (Tat) = 
(jAr_ 2 (c’'i, ■ ■ ■ ,o%, , (Tat) mnltiplied by a fnnction symmetric in (m, v). This is 

trne for 1 < fc < iV. The only remaining case is therefore ai = ■ ■ ■ = Bnt from 
Eq. ( |A.4|) we again have Gn{cfi, ..., ctat-i, <Ji) = Gn- 2 {o' 2 , ■ ■ ■, ctn-i) up to a fnnction 
symmetric in {u,v). We now check explicitly that Gi(cri) and G' 2 (o'i,cr 2 ) are symmetric 
fnnctions in {u,v) for all choices of (cti,(T 2 ). By indnction on N, it follows that Eq. 
( A.4 ) is symmetric in {u, v) for all choices of a^, thereby proving that the inner prodnct 


(^4.38|) is symmetric in (u,v) for all choices of cri,a' G {±s}, as claimed in the text. 


B Relation to the Hofstadter problem 


The Bethe ansatz eqnations (^4.45|) for the case when A^ = 1 have already made a 
surprising appearance [^| in a quantum mechanical context. This is in relation to 
the so-called Hofstadter problem of Bloch electrons on a square lattice in a magnetic 
field. More specifically, in the one-electron problem the Hamiltonian can be expressed 
in terms of the generators of the group of magnetic translations, which in turn can be 
related to 17g(sM(2))-generators at g = (with P,Q mutually prime integers) in 
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the Q = (2s + l)-dimensional representation. The eigenvalue expression and associated 
Bethe ansatz equations were obtained in using the functional Bethe ansatz. In this 
Appendix we recover these equations from our general N results. On the one hand, it 
is a check on our results. On the other, it serves to further elucidate the connection 


found in Ref. [E2 


Dehne then the “Hamiltonian” Ti by 


Comparison with Eq. 


a 

n = 2itTaR^^’^^\u-sX) 

reveals immediately that is simply 
7-f = 2 i — sA), 


(B.l) 


(B.2) 


where the auxiliary transfer matrix on the rhs is for lattice size iV = 1. On the other 
hand, it is clear from Eqs. (|3.1|) and (p.l|) that takes the form 


W = 2ismh(A) (s++ S") 


(B.3) 


in terms of Uq{su{2)) generators. From the eigenvalue expression Eq. ([4.44|) for the 
auxiliary transfer matrix, we hnd the eigenvalue E of to be (after a shift in Uk) 

., sinh[i(M — ua — A)1 
E = 2ismh(M + sA)n -r^TTT-^ 


, ,, sinh[A(M — Un + A)] 

2 i sinh(« - sA) n -tw y - 


(B.4) 


together with the Bethe ansatz equations 


sinh(-Ufc + sA) _ sinh[A(Mfc - uj + A)] 
sinh(Mfc — sA) M sinh[|(Mfc — Uj — A)] 


(B.6) 


Writing q = e^, the Bethe ansatz equations (|B.5|) can be written in the alternative form 


^ 2s 2s ^ 


n 


, fc = 1,..., 2s. 


(B.6) 


g2Ukq-2s _ ^ 1^1 gMfc _ gMjqr’ 

It is clear from Eq. (p.3|) that H is independent of u. Evaluating E at u —>■ oo we 
obtain the expression 

2s 

■ “ (B.7) 


i=i 
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The Hamiltonian Ti has the following interpretation: There exists a g-difference 
operator realization of the (2s + l)-dimensional representation of Uq{su{2)) given by 


= q *T(gz) (B.8) 

= q^^{q~^z) (B.9) 

S+^{z) = -g-2^T(gz)) (B.IO) 

S-^{z) = -z-\q-q-^)-^ (^{q-^z)-'^{qz)) , (B.ll) 

with T(z) belonging to the space of polynomials of degree 2s. In the q —>■ 1 limit, 
this reduces to the familiar differential operator realization of sm(2). In the g-difference 
operator realization, the eigenvalue equation for the Hamiltonian Ti can be written as 

i{zq‘^^ - z-^)'^{q-^z) + i{z-^ - zq-^^)'^{qz) = E'^{z). (B.12) 


Let g = with Q = 2s + 1 and consider hrst the case when P is odd. 


g2« = 


-g ^ and the difference equation ( |B.12|) becomes 

- i(^g“^ + z~^)^{q~^z) + i{z~^ + zq)'<^{qz) = .BT(z), 


We have 
(B.13) 


while the Bethe ansatz equations determining the eigenvalue E become 

Q-l 

Zkq- z 


d + Q 


n 


k= 


(B.14) 


with Zj = e“L 


4(1 + 1 f=i Zk - zjq 
When P is even we have q'^^ = q~^ and we will again arrive at the 


difference equation (|B.13|) if we replace Ti. in Eq. 


hjH = 2isinh(A) - ^+). 


0 1 
1 0 


by 


0 -1 

1 0 


This corresponds, in Eq. (p.l|) , to replacing the twist matrix 

which is also an integrable twist. By Eq. (|4.50|) , we find that the eigenvalue expression 
is i times Eq. (|B.7|) with the same set of Bethe ansatz equations ( p.6|) . The Bethe 
ansatz equations simplify, upon setting e"-’ = izj, to the equations ( p.l4|) without the 
negative sign in front of the product on the rhs. The odd and even P Bethe ansatz 
equations can therefore be unified as 


zj + q 

zlq + 1 


Q-l 

^ qQ TJ ^kQ Zj 

jJi Zk - Zjq'‘ 


with the eigenvalue expression becoming 


/c = 1,..., Q - 1, 


Q-l 


E = -i(g-g ^ 


(B.15) 


(B.16) 


i=i 


The difference equation (p.l3|) is one form of Harper’s equation, which appears in 
the Hofstadter problem with magnetic flux d) = 27 rL per plaquette. Apart from a sign 
factoifl in the Bethe ansatz equations, the solution ( [B.16|) and (|B.15|) for this problem 
was hrst obtained in using the functional Bethe ansatz. 

§We are in agreement with Ref. on the sign discrepancy. 
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